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Abstract 

Quantum Iterated Function System on a complex projective space 
is defined by a family of linear operators on a complex Hilbert space. 
The operators define both the maps and their probabilities by one 
algebraic formula. Examples with conformal maps (relativistic boosts) 
on the Bloch sphere are discussed. 

1 Introduction 

Iterated Function Systems [l] generate fractal sets due to non-commutativity 
of maps. In quantum theory position and momentum operators do not com- 
mute (which leads to the Heisenberg's uncertainty relations), and also differ- 
ent components of spin do not commute. This suggests that fractal patterns 
and chaos may arise as a result of certain quantum measurement processes, 
for instance in a continuous monitoring of several noncommuting observables. 
A typical continuous monitoring takes place, for instant, in a cloud chamber. 
Different regions of the chamber are active in parallel, and they are acti- 
vated sequentially, each at a different time, by a charged quantum particle 
that leaves the track. Moreover, in Heisenberg's picture, position operators 
at different time do not commute. Parallel arrangement of an experimental 
setup is realized by the addition of operators, while serial arrangements lead 
to multiplication of the corresponding operators |2j. While sums of operators 
are commutative, their products, in general, depend on the order of factors. 
Sums of terms usually appear in time evolution generators. Products appear 
in transitions from one quantum state to another, as the result of the wave 
function collapse resulting from measurement events. Repeated application, 
with non-commuting operators, leads to iterated function systems, chaos and 
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fractal attractors on the manifold of quantum states. In the present paper 
we will describe the present status of this new research avenue, and we point 
at some open questions. 



2 The two— sphere S 2 



as the canvas 



The flagship example of a classical "Iterated Function System" (in short: 
IFS) is the Sierpinski fractal. It is generated by an application of 3 x 3 



Figure 1: The classical fractal: Sierpinski Triangle generated by an Iterated 
Function System. 

matrices A[i], i = 1, 2, 3, in a random order, to the vector: 



and ax\ = 0.0, ayi = 0.0, ax 2 = 0.5, ay 2 = 0.0, ax 3 = 0.25, ay 3 = 0.5. (Our 3x 
3 matrices encode affme transformations (of a two-plane) - usually separated 
into a 2 x 2 matrix and a translation vector.) At each step one of the three 
transformations A[i],i = 1,2,3 is selected with probability p[i] = 1/3. After 
each transformation the transformed vector is plotted on the (x, y) plane. 

The important property of the maps A[i] is that they are contractions. 
Sierpinski triangle, as well as another well known example, the fern pQ, live 
on a 2-dimensional plane. Quantum iterated function systems (QIFS) live on 





(1) 



where A[i] is given by 
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complex projective spaces, the simplest one being CP(1) - a 2-dimensional 
sphere S 2 , also known as the Bloch sphere. 

In fact, there are at least two ways in which S* 2 is important in physics. 
First, as the Bloch sphere, it represents pure states of the simplest quantum 
system - spin 1/2. Second, it represents directions in our three-dimensional 
space. The last statement is not relativistically invariant. But there is an- 
other, relativistically invariant interpretation of S 2 , namely as the space of 
directions of light rays. We will start with this second interpretation. 



2.1 S 2 as the projective light cone 

Affine transformations form a natural group of transformations acting on 
the plane. What is the natural group of transformations acting on the two- 
sphere? One would think it is the rotation group 0(3). But rotations are vol- 
ume preserving and they would not mimics contractions. The next candidate 
in line is the Lorentz group 0(3, 1). It is not so well known that the Lorentz 
group acts on a manifold that is diffeomorphic to the sphere S 2 in a natural 
way. One way to see that this is the case is to notice that the Lorentz group is 
the group preserving the space-time metric s 2 = — x^ + x\ + x\ + x\ and thus 
the light cone C = {x = (x°, x 1 , x 2 , x 3 ) : -(x ) 2 + (x 1 ) 2 + (x 2 ) 2 + (x 3 ) 2 = 0}, 
and therefore, because it acts by linear transformations, also the projective 
light cone PO, that is the set of equivalence classes x : x G O, with respect to 
the equivalence relation R C (O\{0}) x (C\{0}) where xRy iff x — Ay, A ^ 0. 
Each equivalence class has a unique representative with x° — 1, so PO can 
be identified with the sphere S 2 = {n G 1R 3 : n 2 = 1}. The Lorentz group 
0(3, 1) consists of 4 x 4 real matrices A = (A^) satisfying A T rjA = r], where 
V = (V(ii>) — diag(— 1, 1, 1, 1) is the diagonal metric matrix. The action 
S 2 9iih A(n) of 0(3, 1) on S 2 is given explicitly by the formula: 



A* + A)n j 
A° + A°.^ ' 



AW 1 = T^-rf-, 0) 



(we use Einstein's summation convention over repeated indices). The group 
0(3, 1) has four connected components. We will need only the connected 
component of the identity S0+(3, 1) consisting of those matrices A in 0(3, 1) 
for which det(A) = 1 and A° > 0. The group SO + (3, 1), though connected, 
is not simply connected. Its simply connected double covering group is the 
group SL(2, C) of 2 x 2 complex matrices of determinant 1. By polar decom- 
position every matrix A G SL(2, C) can be uniquely decomposed A = PU 
into a positive part P and a unitary part U G SU(2). 1 The group SU(2) is 

1 In Relativity the positive matrices represent "Lorentz boosts." 
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the double covering of the rotation group 5*0(3). Nontrivial positive matrices 
in SL(2, C) have two eigenvalues Ai < 1 and A 2 = 1/Ai > 1. It is the positive 
matrices in SX(2,C) that will generate our iterated function systems. 

To describe the 2 : 1 group homomorphism A A(A) from SX(2, C) to 
S0+(3, 1), and also to describe algebraically the action of SL(2, C) on S 2 it 
is convenient to use the Pauli spin matrices a , ax, cr 2 , cr 3 defined by 

1 \ / 1 \ / -i \ ( 1 



= 1 = 1 o i ) 01 = \ i o ) a2 = \i o ; as = v o -i 

The homomorphism SL(2, C) — > ^©-(-(S, 1) is then given by the formula: 

KAf u = \Tr{a»Aa v A*), (4) 

where A* denotes the Hermitian conjugate of A. Every Hermitian 2x2 
matrix X can be uniquely represented as X = x^a^, with x M real. For every 
e G [0, 1], and every unit length vector neS 2 let 

P(n, e ) = l(/ + ea(n)), (5) 

where c(n) = nVi + n 2 o"2 + n 3 cr 3 . It is easy to see that a Hermitian matrix 
X ytz J is an idempotent if and only if it is of the form X = P(n, 1) for 
some n G S* 2 . We will write P(n) = P(n, 1). It is also easy to check that a 
matrix P is positive if and only if it is of the form P = c P(n, e), for some 
c > 0, e G [0, 1], n 6 S 2 . Notice that det(P) = 1 if and only if e < 1 and 
c = 2(1 — e 2 ) -1 / 2 . We will use the matrices P(n, e), with the same e but 
different vectors n to generate IFS-s on S 2 . The formula © describing the 
action of the Lorentz group on S 2 is not the most convenient one for our 
needs. Another way of describing the same action is by noticing that, for 
r G S 2 , the following identity holds 2 

P(n, e)P(r)P(n, e) = A(e, n, r)P(r'), (6) 

where A(e, n, r) > is given by: 

w \ 1 + e 2 + 2e(n • r) 

A(e, n, r) = —± '-, (7) 



while 



„ , (l-e 2 )r + 2e(l + e(n-r))n 



S 2 3 r' = ^ / \ \ " (8) 

1 + e 2 + 2e(n • r) v ' 



2 A more general formula is discussed in Sec. 15.41 Eq. (|23Jl 
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where (n • r) denotes the scalar product 

n ■ r = nin + n 2 r 2 + n 3 r 3 . (9) 

3 The map r i— > r' is the same as the one described in Eq. (J3J), with 
A = A (2P(n,e)/ v / l - e 2 ) . Notice that the dilation coefficient 2/y/l - e 2 
is not important here, because it would cancel out anyway in Eq.©. The 
transformation x^ i— > A.^ u x u implemented by A can be explicitly described by 
the formula known from texts on special relativity: 

x ' = cosh(a)x° + sinh(a)(x • n), (10) 
x' = x — (x • n)n + [sinh(a)x° + cosh(a)(x • n)]n, 

where the "velocity" j3 = tanh(a) = 2e/(l + e 2 ). 4 What is important for 
us, is the fact that the coefficient A(e, n, r) in Eq.© is positive, and thus 
can be interpreted as a (relative) probability associated with the transfor- 
mation r i — > r'. In other words: relative probabilities associated to maps 
implemented by P(n, e) are naturally associated with the maps. It should 
be noticed that positivity of A is guaranteed by the algebraic properties of 
the operators involved. Indeed, because P(n, e) = P(n, e)*, and because 
P(r) = P(r)* = P(r) 2 is an orthogonal projection, the right hand side in 
Eq.© can be represented as A* A, with A = P(r)P(n, e), and is, therefore, 
automatically positive. 



3 Quantum Iterated Function Systems on S 2 . 

Given a sequence i — 1, ... , N of vectors in S 2 we associate with this 
sequence an iterated function system {wi,pi} on S 2 , with place dependent 
probabilities, defined as follows: 

Wt[r) l + e 2 + 2e(n J -r) ' 



A(e, n^r) 
Ef=i A(e, nj, r) 



ftW = ^; V " 7 7 (12) 



3 Because of the property A(e,n, r) > 0, A will be later on - cf. Eq. (fl^|l and Sec. 
15.41 - interpreted as the probability of a jump. This interpretation is natural within the 
quantum measurement theory, but its physical meaning within the framework of light 
directions and Lorentz boosts is not clear. 

4 Notice that for e — > 1, (3 — > 1 - the velocity of light. In this limit the maps r i— > r' 
degenerate to r i— > n and become non-invertible. 
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System {wi,pi} defined by these formulae will be called a Quantum Iterated 



Function Systems or QIFS . 
whenever 



The formula for the probabilities simplifies 



N 

E 



n, 



0. 



(13) 



In the following we will always assume that the vectors n, defining the trans- 
formations Wi add to zero. In this case Pi are given by: 



1 + e 2 + 2e(n, 
iV(l + e 2 ) 



(14) 



In [3] we examined QIFS corresponding to several most symmetric configura- 
tions, where the vectors rij were placed at the vertices of regular polyhedra: 
tetrahedron (4), octahedron (6), cube (8), icosahedron (12), dodecahedron 
(20), double tetrahedron (8), icosidodecahedron (30). In each case numerical 
simulation of the Markov process, starting with a random original point, lead 
to fractal-like patterns on the sphere. For e close to 1 the operators P(n, e) 
are close to projections, therefore the attraction centers are very distinctive. 
For e close to the operators P(n, e) induce maps close to the identity map 
- the patterns are fuzzy. Typical patterns are shown in Fig. 1. 




Figure 2: Quantum Dodecahedron (e = 0.78) and Quantum Octahedron 
(e = 0.58) The darker the place, the smaller probability of it being visited. 



It seems that the fractal dimension depends on the value of e. The Haus- 
dorff dimension of the limit set, for the tetrahedral case, has been numerically 
estimated in Ref. 0] and shown to decrease from 1.44 to 0.49 while e increases 
from 0.75 to 0.95. 



5 A short justification for the term "quantum" will be given in the closing section of 
this paper. 
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4 Transfer (Markov) Operator and Invariant 
Measure 



Stenflo ;5J gives a useful, brief review of the problem of the existence and 
uniqueness of the invariant measures, which is quite useful in our case. We 
will follow the notation and the terminology of j^J. The transfer operator T 
for the system is defined by the formula: 

TV 

(r/)(r) = 5>(r)/K(r)), (15) 

i=l 

where / G C(S 2 ) - the space of all continuous functions on S 2 . By the Riesz 
representation theorem T induces the dual operator T* : \i \— > T*fi on the 
space M(S 2 ) of Borel probability measures on S 2 via the formula: 



/ (T/)(r)d A *(r)= / /(r)d(T»(r). 
Js 2 Js 2 



Since S 2 is compact, there always exists an invariant probability measure 
fj, that is invariant, i.e. T*/x = \i. Numerical simulations of QIFS seem to 
indicate that such a measure is also unique, and that it is concentrated on 
a unique attractor set, though different for different e G (0, 1). As each of 
the normalized operators 2P(n i , e)/yl — e 2 G SL(2, C) has two eigenvalues, 
(1 + e)/(l — e), and (1 — e)/(l + e), - the standard contraction arguments do 
not apply. Nevertheless we have the following theorem: 

Theorem 1. For the Quantum Octahedron the invariant measure is unique 
in the whole range of parameter < e < 1 . 

Proof In Ref. jS] Stenflo states the following classical results, attributed 
to Barnsley et al. Let {(X, d), Pi(x), Wi(x), i G S = {1, 2, . . . , iV}} bean 
IFS with place- dependent probabilities, with all Wi being Lipshitz continuous, 
and all pi being Dini- continuous, and bounded away from 0. Suppose 

/ \ 1 f d(wi(x),Wi(y))\ . . 

sup 2^ Pi (x) log — r < 0. (16) 

x^y i=l \ u>\x, y) J 

Then the generated Markov chain has a unique invariant probability measure. 
The log-average contraction condition ()16j) is somewhat more general than 
the average contraction condition 

Max(w) = supf^ ^f^i < 1. (17) 
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In our case Wi(x) and Pi(x) are analytic, with p*(r) > j^j^s We made 
a numerical estimation of the LHS of the inequality 1)1 7]). with d being the 
natural, rotation invariant, arc distance on S" 2 , for the Quantum Octahedron, 
with < e < 1, and obtained the epsilon-dependence showed in Fig. El, thus 



0.6 - 

H 

0.4 - 

0.2 - 



0_£ 0.4 0.5 0.8 1 



Figure 3: e-dependence of the average contraction parameter Max(w) for 
the quantum Octahedron 



assuring the uniqueness of the invariant measure in this particular case.D 

Let fio be the natural, rotation- invariant, normalized measure on S 2 . 
Then, for any finite n, the measure T* n /i is continuous with respect to 
/io and therefore can be written as 

T*> (r) = /n(r)/i (r). 

The sequence of functions f n ( r ) gives a convenient graphic representation of 
the limit invariant measure. In our case, the functions f n can be computed 
explicitly via the following recurrence formula: 



it 



/n + i(r) = (w-\r)) " f n (wr\r)) 



i=l 



or, explicitly: 



f (r) ^-^ f AKHr)) 

7n+U J iV(l + e 2 ) ^ (1 + e 2 - 2en, ■ r)* UJj 
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where 

! . . (1 — e 2 )r — 2e(l — erij ■ r)ru . , 

1 v 7 1 + e 2 - 2erij • r v ' 

Fig. (J3J) shows a plot of log(/s(r) + 1) for Quantum Octahedron, e = 0.58, 
using the stereographic projection nnz = from S* 2 to the complex 

plane. It should be noticed that via the stereographic projection the maps 
r i— > tWj(r) become fractional, and thus conformal, transformations of the 
complex plane: z \— > = ff^, with a = 1 + enf , 6 = e(n| — in 2 ), c = 

efn, 1 + in 2 ), d = 1 — enf. 




Figure 4: An approximation of the invariant measure: Plot of /s(r) for Quan- 
tum Octahedron (e = 0.58) 



5 Concluding Remarks 

In this section we will place QIFS within a larger field of piecewise-deterministic 
Markov processes and their connection to dissipative dynamics of mixed 
quanto-classical dynamical systems. 

5.1 Classical dynamics 

Usually classical dynamics is described by a 1-parameter group <ftt of diffeo- 
morphisms of a smooth manifold X. In classical mechanics X is a symplectic 
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manifold, the "phase space" of the system, and the flow <p t is generated by 
a Hamiltonian vector field on X. States of the system are simply points of 
X, statistical states are probabilistic measures on X. The set of all statisti- 
cal states is convex, its extremal elements are called pure states. These are 
Dirac measures - concentrated at points of X. The flow <p t on X gives rise to 
a flow on the space of "observables" , that is functions on M, and to a flow 
on the space "statistical states", that is on the space M(X) of probabilistic 
measures on X. If X is discrete, then we can't have a continuous flow on 
X, but we can still have a continuous family of transformations acting on 
observables and on statistical states. 

5.2 Quantum dynamics 

Quantum theory is usually formulated in terms of linear operators acting 
on a separable complex Hilbert space Sj. Observables are represented by 
Hermitian elements of the algebra 21 = L(Sj) of all bounded linear operators 
on ft. Statistical states are positive, normalized, ultra- weakly continuous, 
functionals on 21. They are represented by positive, trace class operators 
p, Tr (p) = 1, with p(A) = Tr(pA), A e 21. Pure states are represented 
by p of the form p = P, where P is an orthogonal projection onto a 1- 
dimensional subspace of f). The space of pure states can be thus identified 
with the space of one-dimensional subspaces of fj. If S) is finite-dimensional, 
f) C n , then the space of pure states is the complex projective plane CP n_1 . 
Quantum dynamics is usually described in terms of a 1-parameter group of 
unitary operators U(t) : t G R. It acts on observables via automorphism 
a t : A ^ U(t)- x AU(t). 

5.3 Mixed quanto-classical dynamics 

We will consider the simple case, where the classical system is finite X = 
{1, . . . , N}. For each a £ X consider the Hilbert space $j a = C Tla and let 
M(n a ) be the algebra of n a x n a complex matrices. 6 The observables of the 
coupled system are now functions a i— > A a e M(n a ) on X with values in 
M{n a ). A pure state of the system is a pair (a,P), where a G {1, . . . , N} and 
P is a Hermitian projection matrix onto a one-dimensional subspace in C ria . 
It is not possible to couple the classical and the quantum degrees of freedom 
via reversible, unitary dynamics. A 1-parameter semi-group of completely 
positive maps of the algebra 21 = ©Q =1 M(n a ) is being used instead. We are 

6 In all examples studied so far the dimensions n a were the same for all a. But such a 
restriction is not necessary. 
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interested in semi-groups with generators of Lindblad's type (also known as 
"dynamical semigroups"), in particular with generators of the form: 

L{A) a = i[H a , A a ] + ^g%A p g Pa - ^(A a A a + A a A a ), (21) 

where g Pa e L(S) a ,Sj/3) and 

A a ^J2^9p a &L(Sj a ). (22) 

We always assume that the diagonal terms vanish: g aa = 0. It has been shown 
in [7] that there is a one-to-one correspondence between semigroups with 
generators of the above type and piecewise-deterministic Markov processes 
on the space of pure states of the system. 



5.4 From dynamical semigroups to QIFS 

Here we are not concerned with the continuous time evolution between jumps, 
so let us extract from the Ref. [7j, and also slightly reformulate, the jump 
process alone. It is determined by the operators g a p alone, and it is an 
iterated function system, with place-dependent probabilities that are also 
determined by g a /3-s. Let (a, P) be pure state, with P being an orthogonal 
projection on a unit length vector ip e Sj a . Observe that for each (3 ^ a we 
have: 

gp a Pg% = \{a,f3-P)Q (23) 

where 

\(a,/3;P)= ||WU 2 >0 (24) 

and, if A(a,/3; P) > 0, then Q is a projection operator on the vector 
gpa^/Wgpai^W in ^3/3- The probabilities p(a,/3; P) are defined as 

p(a, 0; P) = ^ . (25) 

Assume now that all Hilbert spaces S) a = are identical. Assume that 
X = 2 N - the set of N bits, and that g a p = gi ^ when a differs from /3 only 
at one, the z-th bit, otherwise g a p = 0. We will just have a family of operators 
gi and a jump process on pure states P, that is one-dimensional orthogonal 
projections in The maps and their probabilities are determined by: 

gi Pgt = \(i;P)Q, (26) 
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with X(i;P) = ||^|| 2 , Pi (P) = \(i;P)/'E j \U,P), and Q being the or- 
thogonal projection on the subspace spanned by the vector g(ijj. We have an 
iterated function system on the complex projective space CP(n — 1) (equiva- 
lently, on the grassmannian of one- dimensional subspaces of C n ), with place 
dependent probabilities. Both, maps and probabilities, are determined by 
the set of linear operators g i: i £ {1,2,..., N} 



5.5 A short history of QIFS 

The idea of coupling a classical and a quantum system via dynamical semi- 
group has been originally described in jH]. The first model of a QIFS, on 
S 2 , with index a being also continuous a = n, with values in S 2 and, using 
the notation of Sec. g n defined as g n = exp (z7rcr(n)) . These maps were uni- 
tary, thus measure preserving, and did not give rise to a fractal attractor. 
The tetrahedral model was first introduced in ^Uj. It was then examined 
analytically and modelled numerically in a PhD Thesis by G. Jastrzebski 
[3]. The model was further exploited in jllj . where it has been described 
in some details, and where the Lyapunov exponent of the semigroup gener- 
ator has been computed. The term QIFS has been introduced about that 
time on sci.physics.research newsgroup on internet. Recently the term QIFS 
has been adopted in ^2] for another class of maps, namely for maps on the 
space of all statistical states of a quantum system, with arbitrarily assigned 
probabilities. 



5.6 Some open problems 

The question of uniqueness of the invariant measure for a general QIFS is 
an open technical problem. For a particular case of quantum octahedron 
we used numerical estimations of the average contraction parameter. We 
do not have an analytical proof even in this particular case. Then there is 
a "philosophical" question: can the fractal patterns derived from the QIFS 
algorithm be "observed"? Or are they purely mathematical constructs that 
have no relation to the "real world"? The question is not an easy one to 
answer, because QIFS live in the projective Hilbert space of pure states of a 
quantum system, not in "our space". The question of whether they can be 
observed is related to the question: can the wave function be observed? A 
preliminary discussion of this problem has been given in Ref. [Oj. We hope 
to return to this question in our future publications. 
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